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1  Introduction 


This  report  is  concerned  with  optimal  space-time  adaptive  processing  (STAP)  in  environments 
which  cannot  support  fully- adaptive  processing  due  to  either  size,  weight,  and  power  restrictions  or 
sample  support  limitations.  Typical  sensor  platforms  which  exhibit  such  restrictions  include  those 
that  operate  in  airborne,  space-segment,  and  undersea  environments.  The  primary  difficulty  in 
STAP  is  that  the  dimension  of  the  adaptive  weight  vector,  which  equals  the  product  of  the  number 
of  antenna  elements  and  the  number  of  Doppler  bins,  becomes  large.  Both  the  required  sample 
support  for  STAP  detection  and  the  computational  complexity  of  the  algorithms,  which  update 
this  weight  vector,  increase  in  proportion  to  the  dimension  of  the  weight  vector.  When  the  data 
from  auxiliary  range  gates  are  used  for  interference  and  clutter  estimates,  the  required  large  sample 
support  necessitates  the  use  of  range  gates  which  are  relatively  far  from  the  range  gate  of  interest. 
These  distant  range  gates  may  not  preserve  the  homogeneity  of  the  clutter  and  may  not  provide 
good  estimates  of  the  interference  near  the  range  gate  of  interest,  thereby  resulting  in  degraded 
performance. 

Thus,  real-time  STAP  implementations  need  a  lower  rank  processor  to  be  realizable  from  a  com¬ 
putational  viewpoint.  A  dilemma  is  encountered  when  the  optimal  steady-state  solution  requires 
a  larger  number  of  weights  than  the  platform  can  support.  In  this  report  a  general  methodology 
is  developed  which  reduces  the  dimensionahty  of  the  weight  vector  for  target  detection  without 
adversely  affecting  the  steady-state  solution. 

The  optimal  weight  vector  for  adaptive  signal  detection  takes  the  form  of  a  linearly  constrained 
Wiener  filter  [1,  2,  3,  4,  5,  6,  7],  where  the  constraint  is  in  angle  and  Doppler.  It  is  noted  that 
this  linear  constraint  is  not  the  same  as  the  popular  minimum- variance,  distortionless-response 
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(MVDR)  constraint  for  wideband  arrays  [8].  At  first  glance  the  form  of  the  constrained  Wiener 
filter  for  STAP  and  target  detection  does  not  provide  much  insight  to  the  rank-reduction  problem. 
As  a  consequence  there  is  no  clear  methodology  for  finding  the  conditions  needed  to  optimize  a 
partially  adaptive  STAP  system. 

The  most  effective  previously  proposed  reduced-rank  solutions  to  the  STAP  problem  require 
an  eigendecomposition  of  the  array  covariance  matrix  or  its  estimate  [9,  10,  11,  12,  13,  14,  15]. 
These  methods  are  based  on  the  principal  component  approximation  of  a  covariance  matrix  [16, 
17].  One  goal  of  approximating  the  covariance  matrix  with  a  low-rank  estimate  is  to  simplify 
the  weight  update  algorithm.  Notably,  the  weight  vector  obtained  by  reducing  the  rank  of  the 
covariance  matrix  through  the  use  of  the  principal  components  technique  has  full  dimension  [12, 
13,  14,  15].  Finally,  the  principal  component  approach  requires  that  the  dimension  of  the  noise 
subspace  eigenstructure  be  known  precisely.  This  dimension  corresponds  to  the  true  rank  of  the 
interference-plus-clutter  covariance  matrix.  Since  this  rank  is  never  known  a  priori,  it  must  be 
estimated.  The  penalty  for  underestimating  the  rank  of  this  total  noise  contribution  can  be  great; 
the  detection  performance  degrades  rapidly  as  the  rank  of  the  detector  is  decreased  below  the 
dimension  of  the  noise  subspace  eigenstructure. 

In  this  report  a  new  method  is  introduced  to  achieve  a  reduced-rank  array  processor  using 
what  is  called  a  cross-spectral  metric.  This  metric  was  introduced  recently  for  narrowband  MVDR 
processors  and  least  squares  adaptive  filters  [18,  19].  For  these  cases  the  cross-spectral  metric  is 
shown  to  provide  an  upper-bound  on  the  minimum  mean-square  error  performance  obtainable  by 
aU  other  eigen-based  techniques.  In  this  report  the  use  of  the  cross-spectral  metric  is  extended 
to  STAP  for  target  detection.  It  is  shown  that  this  metric  also  maximizes  the  output  signal-to- 
interference  plus  noise  ratio  (SINR)  [6,  20],  thereby  optimizing  the  weight  vector  for  signal  detection 
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Figure  1:  The  Optimal  Joint  Gaussian  Detector. 


under  the  joint  Gaussian  assumption  [3]. 

The  optimal  processor  is  partitioned  by  means  of  a  unitary  operator  which  results  in  a  structure 
called  the  generalized  sidelobe  canceller  (GSC)  [2,  21].  The  GSC  processor  results  in  both  a  lower 
dimensional  weight  vector  and  an  unconstrained  optimal  filter.  The  GSC  is  the  canonical  processor 
for  signal  detection,  and  it  is  by  means  of  this  sidelobe  canceller  interpretation  that  the  process  of 
rank  reduction  is  motivated.  Since  the  cross-spectral  metric  results  in  the  analytical  optimization 
of  the  SINR  as  a  function  of  the  weight  vector  dimension,  a  reduced-rank  weight  vector  is  obtained 
which  best  preserves  the  steady-state  performance. 

The  optimal  processor  for  detection  under  the  joint  Gaussian  assumption  is  interpreted  as 
consisting  of  two  systems  in  cascade,  as  depicted  in  Figure  1.  The  first  system  is  a  noise  field 
estimator  which  utilizes  the  constrained  Wiener  filter.  This  estimation  system  produces  the  best 
estimate  of  the  test  cell  with  which  the  presence  or  absence  of  a  target  is  to  be  determined.  The 
second  system  consists  of  a  hypothesis  testing  mechanism  which  implements  the  detection  criterion. 
The  input  to  the  optimal  processor  is  the  observed  radar  return,  and  the  output  is  a  decision  Di 
with  respect  to  target  presence  {Di  -  1)  or  absence  {Di  =  0).  With  this  interpretation  of  the 
optimal  detector,  the  process  of  partial  adaptivity  belongs  within  the  first  system,  and  results  in 
reduced-rank  estimation. 
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The  space-time  adaptive  processing  framework  is  developed  in  Sect.  2.  In  Sect.  3  the  detection 
criterion  termed  an  AMF  CFAR  test  is  reviewed  with  a  known  covariance  in  terms  of  the  standard 
direct-form  detector.  The  equivalent  fuU-rank  GSC  detector  is  derived  in  Sect.  4.  The  GSC  form 
processor  is  utilized  next  in  Sect.  5  to  derive  a  new  reduced-rank  GSC  AMF  CFAR  detector  which 
uses  the  cross-spectral  metric.  Here  it  is  demonstrated  that  the  cross-spectral  metric  maximizes 
the  output  SINK  as  a  function  of  the  rank  of  the  processor.  The  statistics  of  the  reduced-rank 
AMF  CFAR  test  with  known  and  unknown  covariance  are  derived  in  Sect.  6.  An  example  of 
the  performance  of  partially  adaptive  STAP  for  target  detection  is  presented  in  Sect.  7,  and  the 
conclusions  are  provided  in  Sect.  8. 
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Figure  2:  3- dimensional  CPI  data  cube. 


2  Space-Time  Adaptive  Processing 


Radar  returns  are  collected  in  a  coherent  processing  interval  (CPI),  which  may  be  represented  as 
the  3-D  data  cube  shown  in  Figure  2.  The  data  is  then  processed  at  one  range  of  interest,  which 
corresponds  to  a  slice  of  the  CPI  data  cube.  This  slice  is  a  J  x  K  space-time  snapshot,  denoted  by 
X.  This  space-time  snapshot  can  be  expressed  also  as  the  matrix 


‘'0,0 
I- 1,0 


‘'0,1 


‘'0,/C-l 

h.K-l 


(1) 


where  the  individual  elements  a;j  i;  correspond  to  the  data  from  the  j-th  pulse  repetition  interval 
(PRI)  and  the  A;-th  sensor  element  [3,  22]. 

The  2-D  space-time  data  structure  consists  of  element  space  information  and  PRI  space-Doppler 
information.  If  a  target  is  present  in  the  range  gate  of  interest,  then  the  return  X  is  composed  of 
components  due  to  the  target,  the  interference  sources  or  jammers,  clutter,  and  white  noise: 


X  —  Xt  -b  Xi  -|-  Xc  -f  X^y.  (2) 

If  no  target  is  present,  then  the  snapshot  consists  only  of  interference,  clutter  and  white  noise.  The 
white  noise  in  the  space-time  snapshot  is  the  thermal  noise  present  at  the  sensor  elements,  which 
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is  both  spatially  and  temporally  white.  Each  of  the  other  components  of  the  radar  return  are  now 
addressed. 

The  target  is  assumed  to  approximate  a  point  target,  present  in  only  one  range  gate.  The  J  x  K 
space-time  snapshot  corresponding  to  the  target  is  then  given  by 

Xf  =  (3) 

where  is  the  K  x  1  spatial  steering  vector  in  the  direction  provided  by  the  target  spatial 

frequency  b(wt)  is  the  J  X  1  temporal  steering  vector  at  the  target  Doppler  frequency  Ut,  and 
at  is  the  amphtude  of  the  target.  The  target  data  is  then  stacked  column-wise  to  form  the  AV  x  1 
target  vector 

xt  =  ath{u)t)  ®  (4) 

where  (•)  ®  (•)  represents  the  Kronecker  tensor  product  operator. 

The  interference  environment  consists  of  several  sources  which  are  assumed  temporally  to  be 
white  (barrage  jammers  spread  over  all  Doppler  frequencies  at  a  particular  azimuth).  The  J  X  K 
space-time  snapshot  corresponding  to  the  jammers  is  given  by 

Xi  =  (5) 

where  at  is  the  J  X  I  matrix  of  jammer  amplitudes  for  the  J  PRI’s  of  each  of  the  I  interference 
sources  and  A(t?i)  is  the  K  x  I  matrix  whose  i-th  row  is  the  spatial  steering  vector  for  the  i-th 
interference  source  with  spatial  frequency  The  interference  data  is  stacked  column- wise  to  form 
the  KJxl  interference  vector  x^. 

The  clutter  returns  are  assumed  to  come  from  aU  areas  within  a  common  range  gate  in  a 
direction  towards  which  the  antenna  is  directed.  Consider  returns  from  a  circular  ring  about  the 
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platform.  If  each  spot  within  a  clutter  ring  is  a  point-like  reflector,  its  space-time  snapshot  takes 
the  form, 

Xcp  =  acpb(a;cp)a^(??cp),  (6) 


where  '!?cp  and  Wcp  are  the  spatial  and  Doppler  frequencies  of  the  clutter  return  from  the  p-th  patch. 


and  Ocp  is  the  power  from  the  p-th  clutter  patch.  If  the  clutter  is  assumed  to  be  locally  stationary, 
it  has  a  velocity  relative  to  the  radar  platform  which  is  determined  by  the  azimuth  angle  between 
the  array  and  the  clutter  patch.  Finally,  if  the  clutter  returns  are  assumed  to  come  from  this  clutter 
ring  which  is  divided  into  a  large  number  D  of  evenly  spaced  sectors  and  summed,  the  total  return 
from  the  clutter  has  the  form, 

D 

Xc  =  E  «cpb(/3t?ep)a"(4p),  (7) 

'^Cp=l 

where  (3  is  the  constant  of  linearity  between  the  spatial  and  Doppler  frequencies.  The  spatial 


frequency  is  given  by 


i?cp  =  —cos0csm4>cp, 


with  d,  Ao,  and  6^  being  the  element  spacing,  wavelength,  and  fixed  elevation  angle,  respectively. 
The  clutter  data  is  then  stacked  column- wise  to  form  the  KJ  X  1  vector  Xc. 

The  total  input  noise  vector  n  is  given  by 


n  =  X,  +  Xc  -I-  x^„. 


The  input  noise  covariance  matrix  is  then  defined  to  be 


R  =  E  nn 
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3  The  AMF  CFAR  Test  with  Known  Covariance 


Radar  detection  is  a  binary  hypothesis  problem,  where  hypothesis  Hi  corresponds  to  target  presence 

and  hypothesis  Hq  corresponds  to  target  absence.  Each  of  the  components  of  the  space-time 

snapshot  vector  x  are  assumed  to  be  independent,  complex,  multivariate  Gaussian.  This  snapshot, 

for  each  of  the  two  hypothesis,  is  of  the  form. 

Ho-  x  =  n,  .  . 

Hi:  X  ==  xt  +  n. 

The  KJ  X  1-dimensional  space-time  steering  vector  is  defined  as  follows; 

=  b(a;<)  0  a(7?t),  (12) 

where  h{u>-t)  and  are  introduced  in  (3).  For  convenience  in  the  analysis  to  follow,  the  normal¬ 
ized  steering  vector  in  the  space-time  look-direction  is  defined  to  be 

.  =  ,  '  .  (13) 

The  normalized  target  gain  a  is  then  given  by 

a  =  {■dt,u3t)w{-duut),  (14) 

where  the  original  target  gain  at  is  introduced  in  (3).  The  two  hypothesis  in  (11)  may  now  be 
written  in  the  form 

I;'-  (15) 

:  X  =  as  -f  n,  '  ^ 

where  a  =  |a|  is  a  complex  gain  whose  random  phase  4>  is  uniformly  distributed  between  0 
and  27r.  The  random  vector  x,  when  conditioned  on  is  Gaussian  under  both  hypotheses.  The 


conditional  probability  densities  of  x  are 

/r|Hl,0(x)  = 

A|Ho,.#.(x)  = 


1  (x-as) 

1  ^-x^R-^x 

TT^'^-^iiRir  ’ 
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where  ||  (•)  I1  is  the  determinant  operator.  The  likelihood  ratio  test  then  takes  the  form, 


^  _  f a:\Hi  (x)  _  lo^  /r|ni>(x)#  ^ 

/x|Ho(x) 

where  rj  is  some  threshold.  Using  the  densities  in  (16),  the  test  in  (17)  becomes 


(17) 


A  =  lo  ^2|a|  s^R 


)  77, 

Ho 


(18) 


where  Io(0  is  the  modified  Bessel  function  of  the  first  kind.  The  noise  covariance  matrix  R  is 
nonnegative  definite  and  the  modified  Bessel  function  is  monotonically  increasing  in  its  argument . 
Therefore,  the  test  in  (18)  reduces  to 

,2  Hr 


Ai  = 


s^R-^x 


<  Vu 

Ho 


(19) 


where  the  new  threshold  rji  is  related  to  the  previous  threshold  rj  as  follows: 

‘l-l 


m  = 


2|a| 


(20) 


The  test  in  (19)  was  the  STAP  detection  criterion,  developed  in  the  well  known  papers  by  Brennan 
and  Reed  [3]  and  Reed,  Mallett  and  Brennan  (RMB)  [4].  The  RMB  test  takes  the  form  of  a  matched 
filter  for  incoherent  detection  and  is  a  Bayes  optimal  test  which  is  maximally  invariant  with  respect 
to  the  group  of  phase  shifts.  However,  as  noted  by  Kelly  [23],  no  predetermined  threshold  can  be 
set  to  achieve  a  specified  false  alarm  probability  since  the  detector  is  designed  to  operate  in  an 
unknown  interference  environment. 

The  RMB  test  in  (19)  can  be  modified  to  have  a  constant  frdse  alarm  rate  (CFAR)  property 
via  a  normalization  by  the  output  power.  This  new  test,  derived  independently  in  [24]  and  [25],  is 
termed  an  adaptive  matched  filter  (AMF).  The  AMF  CFAR  test  takes  the  form, 

|s^R“^x 

A2  = 


Hi 
7/0 


(21) 
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and  is  an  optimum  group  invariant  hypothesis  test  with  respect  to  the  change-of-scale  and  the 
initial  target  phase  groups.  The  probability  of  detection  for  the  AMF  CFAR  test  A2  is  given  by 
[24] 

n=0  \  ‘  / 

where  ^  is  the  output  signal-to-interference  plus  noise  ratio,  Q  (•)  is  the  incomplete  gamma  function, 
and  the  threshold  772  is  found  from  the  false  alarm  probability: 

Pfa  =  Qii,V2)  =  (23) 

The  RMB  and  AMF  CFAR  tests  may  be  interpreted  as  functions  of  an  optimal  weight  vector 
WsiNR  [4],  which  is  obtained  by  an  unconstrained  optimization  of  the  SINR  as  follows: 

WsiNR  =  (24) 


where  k  is  some  arbitrary  complex  constant.  The  constrained  optimization  implemented  by  the 
sidelobe  canceller  form,  which  is  introduced  in  the  next  section,  fixes  the  value  of  k  to  provide  a 


normalization  by  the  output  power.  The  weight  vector  in  (24)  is  now  considered  using  a  constrained 
optimization  for  target  detection  under  the  joint  Gaussian  assumption.  The  optimal  constrained 
weight  vector  for  maximizing  the  output  SINR,  while  maintaining  a  normalized  response  in  the 
target  spatio-Doppler  look-direction,  was  originally  given  in  [2]  as  follows: 


WsiNR  — 


R-^s 

s^R-is’ 


(25) 


The  output  y  formed  by  applying  this  latter  weight  vector  to  the  primary  data  is 


y  = 


st^R-i< 


and  the  resulting  output  SINR  is 


^ImR^^SINR 


(27) 
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Figure  3:  The  direct-form  processor. 


where  jap  is  the  output  power  due  to  the  target  signal  under  hypothesis  Hi  and 


s^R“^s 


(28) 


is  the  output  noise  power  of  the  processor.  The  weight  vector  in  (25)  is  the  result  of  a  constrained 
optimization  in  which  a  beam  is  placed  in  that  spatial  and  Doppler  look-direction  which  is  given 
by  the  steering  vector  s.  This  direct-form  processor  is  depicted  in  Figure  3. 

The  RMB  test  in  (19)  may  now  be  expressed  in  the  form 

2  Hi 

I  <  m,  (29) 

which  demonstrates  that  the  RMB  detector  compares  a  threshold  t)i  with  the  ratio  of  the  output 
power  to  the  square  of  the  noise  output  power.  The  normalized  AMF  test  in  (21)  may  similarly  be 
written  in  the  form, 

A2  =  (s^R-^s)  \yf  i  772,  (30) 

Ho 

where  the  threshold  772  is  now  compared  with  the  ratio  of  the  output  power  to  the  noise  output 
power. 
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An  equivalent  new  AMF  CFAR  detector  is  now  developed  which  is  characterized  by  a  sidelobe 
canceller  form.  This  sidelobe  canceller  form  is  used  next  to  reduce  the  dimension  of  the  weight 
vector  WsiNR,  which  maximizes  the  output  SINK  and  the  probabihty  of  detection  under  the  joint 
Gaussian  assumption.  The  goal  of  this  data  reduction  step  is  to  obtain  a  reduced-dimensional  test 
criterion  which  yields  decisions  which  correspond  with  the  fuU- dimensional  AMF  CFAR  test  in 
(21). 
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4  The  GSC  Form  Processor 


The  direct-form  STAP  processor  described  in  Sect.  3  is  now  transformed  into  the  generahzed 
sidelobe  canceller  (GSC)  form.  The  GSC  processor  results  in  an  unconstrained  weight  vector 
and  reformulates  the  STAP  detector  structure  into  the  form  of  a  standard  Wiener  filter.  This  GSC 
form  of  the  detector  provides  a  deeper  insight  to  the  rank-reduction  problem  than  appears  to  be 
possible  otherwise. 

Consider  now  using  a  unitary  KJ  X  KJ  matrix  operator  T  to  transform  the  data  prior  to 
adaptive  detection.  The  unitary  nature  of  this  matrix  conserves  power,  so  that  the  resulting  output 
SINR  is  seen  to  be  identical  to  the  original  direct-form  detection  processor.  The  structure  of  this 
operator  is  partitioned  as  follows; 


T  = 


(31) 


where  the  KJ  x  1-dimensional  conventional  space-time  beamformer  s  is  defined  in  (13)  and  B  is 
the  fuU-row  rank  N  x  KJ  signal  blocking  matrix  which  maps  x  onto  the  null-space  of  s,  where 
N  =  (KJ  -  1).  Hence, 


Bs  =  0, 


(32) 


so  that  the  matrix  B  effectively  blocks  any  signal  coming  from  the  spatio-Doppler  look-direction. 
Any  full  row-rank  matrix  B  that  satisfies  (32)  and  results  in  an  invertible  T  is  a  valid  signal 
blocking  matrix.  The  Gram-Schmidt  algorithm  may  be  applied  then  to  this  matrix  to  generate 
an  orthonormal  B  and  a  unitary  T.  Two  algorithms  for  directly  finding  an  orthonormal  B,  using 
the  singular  value  decomposition  (SVD)  and  the  QR  decomposition,  respectively,  are  described  in 
Appendix  A. 

The  transformation  of  the  radar  return  x  by  the  operator  T  in  (31)  yields  a  vector  x  which  has 
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the  form 


X  =  Tx  = 


s^x 

d 

Bx 

b 

(33) 


where  the  scalar-valued  beamformed  output  is  denoted  by  d.  Here  also  the  iV-dimensional  vector 
b  is  termed  the  noise  subspace  data  vector.  It  is  seen  now  that  the  transformed  data  vector  x  has 
an  associated  covariance  matrix  which  takes  the  form: 


Rj  =  TRT^ 


^bd 


m  Rfc 

The  N  X  N  noise  subspace  covariance  matrix  Rfc  is  expressed  by 


(34) 


Ri,  =  E[bb^]  =  BRB^.  (35) 

The  iV  X  1  cross-correlation  vector  between  the  noise  subspace  data  vector  and  the  beamformer 
output  is  given  by 

rid  =  E[bd*]  =  BRs,  (36) 

where  *  represents  the  complex  conjugate  operator.  The  scalar  in  (34)  is  computed  to  be 

aj  =  s^Rs,  (37) 


and  represents  the  variance  of  the  conventional  beamformer  output. 


Next  let  T  operate  on  the  steering  vector  s.  This  operation  yields  the  unit  transformed,  steering 


vector  ei ,  given  by 


ei  =  Ts  = 


1 

0 


0 


Then  the  optimal  weight  vector  in  these  transformed  coordinates  is  given  by 


^gsc 


efRjei 


1 

— w 


7 


(38) 


(39) 
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Figure  4:  The  full-rank  GSC  processor. 


where  is  the  covariance  matrix  in  (34).  The  partitioning  of  the  matrix  operator  T  leads  naturally 
to  the  form  of  the  GSC  array  processor  depicted  in  Figure  4  [2,  7,  21].  This  processor  results  in  a 
fixed  weight  of  unity  for  the  upper  branch  and  an  adaptive  weight  vector  w  of  dimension  N  =  KJ-l 
in  the  lower  branch.  The  vector  w  in  (39)  is  provided  by  the  Wiener  solution  corresponding  to  the 
filter  depicted  in  Figure  4; 

w  =  (40) 

where  Rj,  and  vu  are  computed  in  (35)  and  (36),  respectively.  The  GSC  form  processor  implements 
the  FT J-dimensional  weight  vector  in  (39)  using  the  partitioning  defined  in  (33).  The  steady-state 
performance  of  the  GSC  and  the  direct-form  processor  are  identical,  but  the  adaptive  weight  vector 
w  in  the  GSC  is  of  a  lower  dimension.  Hence,  the  computational  requirements  for  updating  this 
weight  vector  are  reduced  and  the  GSC  form  therefore  can  be  considered  canonical. 
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A  normal  component  covariance  matrix  Rp,  cross-correlation  vector  Vpd,  and  Wiener  filter  are 
defined  now  as  follows: 


Rp 

=  E[pp^] 

=  U^RfcU  =  A, 

(46) 

^pd 

=  E[pd*]: 

=  U^rtrf, 

(47) 

Wn 

=  Rp^rpd 

=  U^w. 

(48) 
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Figure  5:  The  full-rank  GSC  processor  in  principal  coordinates. 


The  GSC  in  these  normal  coordinates,  depicted  in  Figure  5,  is  equivalent  to  the  GSC  in  Figure  4 
in  terms  of  its  steady-state  characteristics.  The  array  output  of  the  GSC  in  normal  coordinates  is 
given  by 

y  =  X  =  (s^  -  x.  (49) 

Note  that  the  output  noise  power, 


P  =  (^d-  rprfRp  =  <^d-  ^^d, 


(50) 


and  that  the  SINK, 


a 


a 


^d  -  ^  rpd  aj  -  ^  nd  ’ 

are  conserved  by  any  unitary  transformation,  including  that  realized  by  the  operator  U. 


(51) 
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5  Partially  Adaptive  STAP  and  the  Cross-Spectral  Metric 


We  now  derive  the  cross-spectral  metric  for  optimal  subspace  selection  and  rank  reduction  in 
partially  adaptive  STAP.  The  problem  of  reducing  the  degrees  of  freedom  for  an  array  processor 
involves  selecting  a  subset  or  some  combination  of  the  elements  to  be  adaptively  weighted.  For 
notational  purposes,  let  the  space  spanned  by  the  columns  of  the  fully  adaptive  array  covariance 
matrix  be  denoted  by  ,  implying  that  the  observation  covariance  matrices  Rj,  and  Rp  are  of 
dimension  iV  X  iV  and  the  vectors  rpd,  w,  and  are  N  X  1-dimensional  vectors.  The  partially 
adaptive  GSC  shown  in  Figure  6  utilizes  an  iV  X  Af  transformation  operator  U,  in  place  of  U  in 
Figure  5,  to  form  the  M-dimensional  reduced-rank  observation  data  vector, 

z  =  U^h,  (52) 

where  M  <  N.  The  associated  M  x  M  reduced-rank  covariance  matrix  is  given  by 

R,  =  U^RbU  =  Am,  (53) 

where  Am  is  the  diagonal  matrix  composed  of  the  M  eigenvalues  corresponding  to  the  eigenvectors 
to  be  selected  which  form  U.  The  cross- correlation  between  the  process  z  and  the  desired  signal  d 
is  given  by 

T,d  =  E[zd*]  =  U^Tbd.  (54) 

The  data  vector  z  is  then  processed  by  the  reduced-rank  weight  vector  Wm,  which  is  of  dimension 
M  X  1.  From  Wiener  filter  theory  this  weight  vector  is  expressed  finally  by 

Wm  =  R^^r^jd  =  A^^U^Tbd-  (55) 

The  most  popular  technique  for  subspace  selection  is  based  on  the  principal  components  method 
[9,  11,  12,  13,  14].  This  method  determines  the  singular  value  decomposition  of  the  N  X  N- 
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dimensional  covariance  matrix  Rj  and  selects  the  M  largest  eigenvectors  (those  corresponding  to 
the  largest  eigenvalues)  to  form  the  M-dimensional  eigen-subspace  ^  C  in  which  the  adaptive 
processor  operates.  However,  this  technique  does  not  directly  consider  the  maximum  output  SINK 
performance  measure  to  maximize  the  probability  of  detection,  which  is  a  function  of  not  only 
the  space  spanned  by  the  noise  covariance  matrix  R;,  but  also  of  the  cross-correlation  between  the 
desired  signal  d  and  the  noise  process  b. 

The  derivation  in  [18]  is  now  used  to  reduce  the  rank  of  the  Wiener  filter  in  a  manner  which 
maximizes  the  output  SINK.  To  accompUsh  this  a  rank  reducing  N  x  M  transformation  matrix  U, 
which  is  composed  of  some  M  columns  from  U  ,  needs  to  be  selected  (see  Figure  6).  The  operator 
14  is  constrained  therefore  to  be  a  subset  of  M  of  the  N  possible  eigenvectors  of  R(,.  This  particular 
constraint  allows  a  direct  comparison  with  the  principal  component  technique,  which  chooses  the 
rank  reducing  transform  to  be  composed  of  those  M  eigenvectors  which  correspond  to  the  largest  M 
eigenvalues.  Thus,  the  particular  problem  at  hand  is  to  choose  the  subspace  spanned  by  a  set  of  M 
eigenvectors  out  of  the  N  available  such  that  the  resulting  Af-dimensional  Wiener  filter  yields  the 
largest  output  SINR  out  of  all  (j^^)  possible  combinations  of  eigenvectors.  The  subspace  spanned 
by  the  columns  of  the  optimal  reduced-rank  covariance  matrix  which  achieves  this  maximization 
of  the  output  SINR  as  a  function  of  filter  rank  is  denoted  by  1?  C  . 

Now  denote  the  reduced-rank  processor  output  by  t/r-  This  is  illustrated  in  Figure  6.  A  study  of 
Figures  5  and  6  suggest  that  the  reduced-rank  processor  output  may  be  expressed  by  the  relation. 


Vr  = 


1, 


y 

b 


(56) 


Denote  the  weight  error  vector  between  the  fuU-rank  weight  vector  and  its  reduced-rank  version  by 


e  =  -  Uwm- 


(57) 
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corresponds  with  the  selection  of  the  M  columns  of  U  that  were  retained  to  form  the  rank  reducing 
operator  U.  Now,  one  wants  to  minimize  the  scalar  term  e^'R-bS  which  appears  in  (58)  and  (59) 
as  follows: 

min  [e^R6e  =  min  (R^ Vjv  -  JR^ Vm)]  ■  (61) 

Thus  evidently  the  best  solution  for  (61)  is  to  choose  that  set  of  M  rows  of  U  for  the  M  x  N  matrix 
U  such  that  JR^^Wm  is  the  best  low  rank  approximation  to  the  vector  Rp^^Wjv. 

The  Wiener-Hopf  relationship  for  the  full-rank  case  is  given  by 
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for  i  =  1, 2, . .  .iV.  Note  that  the  metric  in  (65)  measures  the  cross-spectral  energy  projected  along 
the  j-th  eigenvector.  With  this  selection,  the  columns  of  the  reduced-rank  covariance  matrix  R,;;; 
span  the  M-dimensional  cross-spectral  subspace  C  to  provide  the  largest  output  SINK  of 
any  M-dimensional  subspace  which  is  spanned  by  M  of  the  N  columns  of  U. 

Clearly,  the  subspace  H  spanned  by  the  columns  of  eigenvectors  corresponding  with  the  M 
largest  values  of  the  cross-spectral  metric  is  not  the  same  as  the  subspace  ^  spanned  by  the 
eigenvectors  corresponding  with  the  M  largest  eigenvalues  for  aU  values  of  M.  This  means  that  the 
Wiener  filter  in  the  cross-spectral  subspace  £2  yields  an  output  SINR  which  is  always  larger  than  or 
equal  to  that  provided  by  the  Wiener  filter  in  the  subspace  'i' .  We  conclude  that  the  reduced-rank 
Wiener  filter  found  via  the  cross-spectral  metric  maximizes  the  output  SINR  as  a  function  of  the 
rank  of  the  filter  for  the  eigenvector  basis.  Thus  the  cross-spectral  metric  results  in  an  upper-bound 
on  the  performance  of  eigen-based  rank  reduction  techniques. 

To  demonstrate  that  the  cross-spectral  metric  is  optimal  for  each  rank  M  <  iV,  we  consider 
the  decomposition  of  the  SINR  performed  by  the  full-rank  matrix  of  eigenvectors  U.  The  output 
SINR  of  the  full-rank  processor,  from  (46)  and  (51),  can  be  expressed  as 


Vfvu 


Finally,  a  comparison  of  the  decision  rule  in  (65)  with  the  expressions  in  (66)  and  (67)  demonstrate 
that  the  cross-spectral  metric  maximizes  the  output  SINR  as  a  function  of  the  rank  of  the  weight 


vector. 
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6  The  Reduced- Rank  AMF  CFAR  Test  with  Known  and  Un¬ 
known  Convariance 


The  AMF  CFAR  test  in  (21)  is  defined  in  terms  of  the  output  signal  y  and  the  output  noise  power 
P  of  the  direct  form  processor.  The  GSC  is  then  derived  as  the  canonical  form  of  this  optimal 
array  processor  in  Sect.  4,  where  it  is  also  demonstrated  that  the  GSC  yields  an  output  signal  and 
SINR  which  are  equivalent  with  that  provided  by  the  direct-form  processor. 

An  optimal  reduced-rank  STAP  design  technique  is  derived  in  Sect.  5,  where  the  GSC  form 
processor  was  utiUzed.  It  is  also  shown  in  Sect.  5  that  the  target  signal  power  is  not  affected 
by  rank  reduction.  Thus,  as  mentioned  in  Sect.  1,  the  process  of  rank  reduction  for  partially 
adaptive  STAP  belongs  within  the  noise  field  estimation  system  shown  in  Figure  1.  The  cross- 
spectral  metric  provides  an  upper-bound  on  the  output  SINR  achievable  by  reduced-rank  detectors 
using  an  eigenvector  basis.  This  fact  explains  why  the  cross-spectral  metric  provides  a  capability 
to  achieve  reduction  in  the  rank  of  the  weight  vector  below  the  dimension  of  the  noise  subspace 
eigenstructure  without  a  significant  loss  in  performance.  The  cross-spectral  metric  given  in  (65) 
provides  the  best  set  of  the  possible  (^)  eigenvectors  of  the  noise  subspace  covariance  matrix  to 
select  for  rank  reduction.  These  optimal  eigenvectors  are  then  used  by  the  hypothesis  testing 
mechanism  to  determine  target  presence  or  absence.  A  new  optimal  reduced-rank  AMF  CFAR  test 
is  now  formally  defined  and  analyzed. 


For  convenience,  we  now  rewrite  the  AMF  CFAR  test  in  (30)  as  follows: 


s^R 

s^R“^s 


IM!  >■„ 

p  <  ^2- 

^  Ho 


The  reduced-rank  array  output  yr,  derived  in  (56)  and  depicted  in  Figure  6,  is  written  now  in  the 
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form, 


yr  =  (s^  -  X  =  =  d-  ^z. 

Also  the  reduced-rank  output  noise  power  in  (58)  may  be  expressed  by 

Pr^(^d- 


(69) 


(70) 


where  the  covariance  matrix  and  the  cross-correlation  vector  r^d  were  defined  in  (53)  and  (54), 
respectively.  In  terms  of  the  expressions  in  (69)  and  (70),  a  new  optimal  reduced-rank  version  of 
the  AMF  CFAR  test  is  obtained  as  follows: 


h.2r 


<  >?2, 

Ho 


(71) 


where  the  reduced-rank  SINK  is  given  by 


crj  -  r^dRz 

The  probability  of  detection  of  this  new  reduced-rank  AMF  test  in  (71)  is  modified  to 

and  the  probability  of  false  alarm  retains  the  form. 


(72) 


(73) 


PFA  =  Q{hm)  =  e-^^- 


(74) 


Note  that  the  expressions  in  (69),  (70),  (71),  (72),  (73)  and  (74)  are  valid  for  any  unitary  operator 
used  for  rank  reduction  with  a  subspace  selection  rule  which  chooses  M  of  the  A  columns  of  the 
unitary  operator.  This  includes,  but  is  not  limited  to  the  reduced-rank  principal  component  and 
cross-spectral  AMF  CFAR  tests,  which  use  the  unitary  matrix  composed  of  the  eigenvectors  of  Rj. 

The  reduced-rank  AMF  CFAR  test  is  now  examined  for  the  case  of  unknown  covariance.  The 
standard  approach  for  radar  detection  when  the  noise  covariance  is  unknown  is  to  obtain  an  estimate 
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from  some  L  >  KJ  range  gates  in  the  neighborhood  of  the  range  gate  currently  being  tested  for 
signal  presence.  The  data  in  the  range  gate  of  interest  is  termed  the  primary  data.  The  L  nearby 
range  gates,  which  yield  what  is  termed  the  secondary  or  auxiliary  data,  are  assumed  to  be  target- 
free.  It  is  also  possible  to  assume  that  only  multiple  snapshots  of  the  same  range  gate  are  available. 
The  presentation  of  the  material  in  this  section  follows  similar  derivations  given  previously  for 
fuU-rank  detectors  [26,  27,  24,  25,  28,  29,  30]. 

It  is  well  known  that  the  maximum  likelihood  estimate  of  the  covariance  matrix  R  is  provided 
by  the  sample  covariance  matrix  R,  i.e. 


where  rij  is  the  KJ  x  I  vector  of  samples  from  the  i-th.  nearby  range  gate.  Define  an  auxiliary 
KJ  X  L  data  matrix  A  by 

A  =  [ni,  n2,  •••  ni,],  (76) 

so  that  the  sample  covariance  matrix  R  may  be  written  in  the  form, 

R  =  yAA^.  (77) 

Then  the  AMF  CFAR  test  is  expressed  in  terms  of  the  sample  covariance  matrix  as  follows: 

The  data  matrix  A  is  expressed  now  in  the  form 

A  ^  (79) 

where  X  is  a.  KJ  X  L  matrix  of  zero-mean,  independent  and  identically  distributed,  complex 
Gaussian  variates  with  unit  variance.  Then  an  equivalent  test  to  the  sampled  data  AMF  CFAR 


L  s"  (aA^)  ^  X  Hi 

- - - - <  ^2 

s^^(aA^)“%  Ho 
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test  in  (78)  is  provided  by 


Ao 


M 


-1 


X 


-1 


Hi 

<  m- 

Ho 


(80) 


The  methodology  used  in  the  previous  section  for  transforming  the  observed  data  into  the  GSC 
form  is  now  applied  to  the  auxiliary  data.  The  matrix  T  defined  in  (31)  operates  on  the  data 
matrix  X  defined  in  (79).  The  transformed  data  matrix  X  has  the  form, 


X  =  TX  = 


V 

B 


(81) 


where  Tf  is  the  1  x  T  beamformed  data  vector  and  B  is  the  N  X  L  noise  subspace  data  matrix 
where  N  =  KJ  —  1. 

The  rank  reducing  N  x  M  transform  U,  defined  to  be  composed  of  those  M  eigenvectors  of  R;, 
which  maximize  the  cross-spectral  metric  given  in  (65),  is  applied  next  to  form  a  (M  -f  1)  x  (M  -f- 1) 


transformed  reduced-rank  sample  covariance  matrix  as  follows: 


1  0^^ 

0 


TXX^T^ 


'  1 

0^  ■ 

H 

T>V^ 

T>Z^  ' 

0 

ZV^ 

ZZ^ 

Rz 

(82) 


where  the  reduced-rank  M  x  T- dimensional  auxiliary  data  vector  Z  is  given  by 


Z  =  U^B. 


(83) 


The  optimal  (M  -f  1)  X  1-dimensional  reduced-rank  weight  vector  in  these  transformed  coordinates 
is  given  by 

(84) 

where  the  M  x  1-dimensional  vector  is  provided  by  the  Wiener  solution  corresponding  to  the 
filter  depicted  in  Figure  4,  namely 


w,  =  [ZZ^)  '  ZV^  = 


(85) 
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The  resulting  output  power  is  given  by 


Pr  =  al- 


(86) 


The  primary  data  is  also  transformed  by  R  and  T  to  obtain 


x„  =  TR  = 


d 

b 


(87) 


and  the  reduced-rank  M  X  1-dimensional  primary  data  vector  is  then  given  by 


Zp  =  U^h. 


(88) 


The  reduced-rank  AMF  CFAR  test  with  an  unknown  covariance  corresponding  to  the  test  in  (78) 
is  expressed  now  in  the  form, 

Id  -  wifzj  \y  |2  Hi 

?  »2,.  (89) 


^2sr  = 


rzd  Pr  Ho 

where  the  threshold  ri2s  in  (89)  is  related  to  the  original  threshold  772  in  (78)  by 


m 

L' 


(90) 


Next  the  denominator  of  the  reduced-rank  sampled  data  AMF  CFAR  test  in  (89)  is  rewritten 
in  the  form, 

Pr  =  V  (1- =  VCV^.  (91) 

The  vector  X>  is  a  1  X  T  zero-mean,  white  Gaussian  random  vector.  Therefore  the  distribution  of 
V  can  be  denoted  by 

V  ^  iv,(0,  l),  (92) 

where  Ai(0, 1)  is  a  joint  Gaussian  density  of  zero  mean  and  covariance  I,  the  LxL  identity  matrix. 
Also  it  is  easily  shown  by  (83)  that  the  matrix  I  —  C,  given  by 


(I  -  c)  =  z^k-^z 


(93) 
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is  a  projection  matrix  of  rank  M  which  projects  vectors  onto  the  reduced-rank  auxiliary  data  noise 
subspace.  The  matrix  C  is  also  an  idempotent  projection  matrix  of  rank  L—M,  and  the  distribution 
of  Pr  is  therefore  central  with  L  —  M  degrees  of  freedom  [31,  Theorem  1.4.2,  p.27].  Thus  the 
distribution  of  Pr  is  expressed  by 

Pr  ^  xLm-  (M) 

The  output  of  the  reduced-rank  GSC  form  array  variate  may  be  expressed  next  by 

2/r  =  d- (95) 


c 


d  ^  pir- 


(100) 
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The  parameter  7  is  the  ratio  of  two  independent  variables  where  the  numerator  is  noncentral 
X^  distributed.  Hence,  under  hypothesis  7  has  a  noncentral  F  distribution  [31,  Theorem  1.3.6, 


p.24]  given  by 


/y|p,Hi(7)  = 


/3(1,X-M)(1  +  7)^'-^+i 


iTi  [  L  -  M  +  1,1, 


6  7 
1  +  7 


(101) 


Therefore,  the  conditional  distribution  of  the  test  A2s^  given  Zp  and  Z  under  hypothesis  Hi,  may 
be  expressed  by 


e  ^  p  f  ^  Q  P 

=  13(1,  L-  M)  (1  +  1^1  +  + 


(102) 


The  unconditional  probability  density  function  of  the  reduced-rank  AMF  CFAR  test  is  then  com¬ 
puted  as  follows: 

/A2j^|ni(9)  =  /  /A23^|p,ni(9) /p|ni(p)  (103) 

J  0 

which  is  expressible  in  the  form. 


/a2„ jHi(g)  - 


/?(l,T-M)/3(M,X-M  +  l)io  (1  +  gp)^-^+i 


I 


1+1  T  —  M  +  1, 1, 


1  +  «'P 
(104) 


where  g  >  0.  Finally  by  (104)  for  a  fixed  threshold  t],  the  probabihty  of  detection  is  computed  to 
be 

roo 

PDiv)=  ijij(g)dg,  (105) 

Jt)  ^ 

which  when  evaluated  has  the  form 


Pd{v)  =  [ 

Jo 


(1  -  pY 


I3{1,L-M)I3{M,L-M  +  1) 


/  (1  -  ^Fi  {L-M +  1,1, 6q)  dq  dp. 

i  +  VP 


(106) 


dp. 
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The  distribution  of  the  reduced-rank  AMF  CFAR  test  under  hypothesis  Hq  requires  that  the 
variable  6  be  set  to  zero,  and  (104)  reduces  to 


/A2.Jno(^)  - 


-1-1)70  (l  +  gp)^-^+i 


I  „  ,  ■  I-mI,  ■  ip.  (m 

Jo 


which,  for  the  selected  threshold  ?y,  yields  a  false  alarm  probability  given  by 


roo  \  roo  t 

Pfa{v)  =  IHo (q)  dq  =  +  1)  X  L 


oo  rl  (1  _  p) 


M-l 


(1  -I-  qp)^-^+^ 


dpdq. 

(108) 


Since  the  PpA  is  only  a  function  of  the  threshold  rj  and  the  integers  L  and  M,  aU  of  which  are  held 
fixed,  the  test  A2s,,  has  a  constant  false  alarm  rate. 
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7  An  Example  of  the  Performance  of  Partially  Adaptive  Radar 


A  simulation  is  now  considered  which  utilizes  the  parameters  of  the  ARPA  Mountain  Top  radar 
in  order  to  compare  the  performance  of  the  reduced-rank  detectors  to  that  of  the  full-rank,  joint 
domain  optimal  detector.  The  AMF  CFAR  test  with  known  covariance  is  utilized  exclusively  for 
this  example. 

The  Mountain  Top  radar  employs  the  Radar  Surveillance  Technology  Experimental  Radar 
(RSTER)  and  the  Inverse  Displaced  Phase  Center  Array  (IDPCA),  co- located  at  the  same  site 
[33].  It  is  assumed  that  the  radar  is  in  the  RSTER-90  configuration  and  the  receive-only  mode. 
The  transmit  frequency  is  450  MHz.  This  radar  consists  of  K=14  elements  and  J=16  pulses  in 
the  PRI.  The  elevation  angle  6  is  fixed  (pre-beamformed),  and  the  azimuth  angle  (/>  is  the  only  free 
parameter.  The  dimension  of  the  adaptive  processor  is  KJ  =  224.  The  GSC  processor  beamforms 
the  target  signal  in  one  branch  and,  using  a  bloching  matrix  B,  blocks  the  target  signal  from  the 
other  branch.  The  noise  subspace  data  vector  is  then  of  dimension  A  =  A  J  —  1  =  223.  The 
corresponding  full-rank  GSC  weight  vector  is  also  of  dimension  223  x  1. 

For  the  purposes  of  this  analysis,  three  barrage  jammers  and  land  clutter  compose  the  inter¬ 
ference  enviroment  and  one  target  signal  is  present.  It  is  assumed  that  the  radar  platform  is  at  a 
height  of  500  meters,  and  that  the  platform  velocity  is  500  meters  per  second.  The  target,  which 
is  also  at  a  height  of  500  meters,  has  a  range  of  22  kilometers  and  an  azimuth  of  -30°.  The  target 
velocity  is  250  meters  per  second.  This  target  has  a  signal-to-noise  ratio  (SNR)  of  15  dB,  represent¬ 
ing  a  small,  non-fluctuating,  constant  radar  cross-section  target.  The  three  jammers  have  azimuth 
angles  of  —60°,  30°  and  60°,  with  jammer-to-noise  ratios  of  30  dB,  40  dB,  and  30  dB,  respectively. 
The  clutter-to- noise  ratio  is  set  to  be  30  dB. 
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The  cross-spectral  metric  is  now  applied  to  the  GSC  and  the  detection  performance  is  compared 
to  the  principal  component  techniques  such  as  PCI  [13,  14],  which  has  also  been  termed  an  eigen- 
canceler  [12,  15].  The  performance  of  the  fully  adaptive  GSC,  the  partially  adaptive  eigen-subspace 
GSC,  and  the  partially  adaptive  cross-spectral  subspace  GSC  processors  are  evaluated. 

To  obtain  the  performance  of  the  subspace  selection  techniques,  let  the  dimension  of  the  adaptive 
processor  be  reduced  from  N  =  223  weights  to  M  =  50.  The  dimension  of  the  noise  subspace 
eigenstructure  is  calculated  to  be  approximately  70.  This  is  the  lower  bound  for  the  standard  rank 
reduction  techniques  [9,  10,  11].  The  power  spectrum  of  the  data,  averaged  over  500  snapshots, 
is  presented  in  Figure  7,  from  which  it  is  evident  that  the  target  signal  is  not  discernible.  Figure 
8  presents  the  gain  of  the  full-rank  spatio-temporal  Wiener  filter,  from  which  it  is  seen  that  the 
optimal  space-time  weight  vector  of  dimension  223  is  capable  of  attenuating  all  jammers  and  clutter 
while  passing  the  target  signed.  The  full-rank  Wiener  filter  SINR  is  38.28  dB.  Figure  9  depicts  the 
spatio-temporal  gain  of  the  rank  50  eigen- subspace  Wiener  filter.  The  clutter  and  most  of  the  power 
from  the  three  jammers  is  passed  with  relatively  high  gain;  there  is  little  attenuation  of  the  noise. 

The  output  SINR  for  the  eigen-subspace  Wiener  filter  is  25.98  dB,  reflecting  a  loss  of  12.3 
dB.  The  spatio-temporal  gain  of  the  rank  50  cross-spectral  subspace  Wiener  filter  is  depicted  in 
Figure  10.  It  can  be  seen  that  by  simply  selecting  a  different  set  of  50  eigenvectors,  the  resulting 
array  response  closely  approximates  the  fuU-rank  optimal  STAP  matched  filter  and  the  clutter  and 
jammers  are  attenuated.  The  resulting  output  SINR  is  28.03  dB,  reflecting  a  loss  of  only  0.25  dB 
in  reducing  the  rank  from  223  to  50. 
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Normalized  Doppler 
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Figure  9:  The  reduced-rank  eigen-subspace  Wiener  filter. 
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Normalized  Doppler 
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Figure  10:  The  reduced-rank  cross-spectral  subspace  Wiener  filter. 
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Weight  vector  rank 

Figure  11:  The  output  SINR  of  the  reduced-rank  processors  as  a  function  of  the  rank  of  the  data. 
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The  performance  of  these  two  reduced-rank  processors  may  also  be  evaluated  by  plotting  the 
output  SINR  as  a  function  of  Wiener  filter  rank.  The  SINK  for  the  eigen-subspace  and  cross- 
spectral  subspace  reduced-rank  processors  are  compared  to  the  full-rank  SINR  in  Figure  11  as  the 
rank  increases  from  1  to  223.  It  can  be  seen  from  this  figure  that  the  eigen-subspace  processor 
converges  to  the  optimal  SINR  in  an  abrupt  manner,  thereby  obtaining  the  optimal  value  only 
once  the  Wiener  filter  rank  reaches  the  dimension  of  the  noise  subspace  eigenstructure.  This  is 
characteristic  of  the  eigen-subspace  techniques  in  nearly  aU  of  its  applications  since  the  principal 
components  method  does  not  yield  an  optimal  filter  as  a  function  of  rank.  By  contrast  the  cross- 
spectral  metric  depicts  an  exponential-like  convergence  to  the  optimal  SINR. 

The  above  simulation  shows  that  the  reduced-rank  Wiener  filter  in  the  cross-spectral  subspace 
outperforms  the  Wiener  filter  in  the  eigen-subspace  for  aU  filter  ranks  less  than  the  dimension 
of  the  noise  subspace  eigenstructure.  It  is  verified  in  Figure  11  that  the  cross-spectral  metric 
provides  an  upper-bound  to  the  SINR  performance  achievable  by  a  reduction  of  the  number  of 
eigenvectors  of  the  array  covariance  matrix.  FinaRy,  note  for  the  Mountain  Top  radar  example  that 
the  reduction  of  the  rank  to  50  weights  can  be  interpreted  as  an  error  in  estimating  the  dimension  of 
the  noise  subspace  eigenstructure.  It  is  clear  from  Figure  11  that  the  penalty  in  SINR  performance 
for  underestimating  the  rank  of  the  noise  subspace  eigenstructure  is  greatly  reduced  in  the  cross- 
spectral  approach  to  rank  reduction.  This  robustness  property  is  paramount  for  reduced-rank 
detection  in  unknown  noise  environments. 

Next  the  probability  of  false  alarm  for  each  test  is  set  to  PpA  =  10“^.  The  detection  probabilities 
are  calculated  for  the  full-rank  test  (223  adaptive  weights)  and  for  the  two  reduced-rank  tests  (50 
adaptive  weights)  in  the  eigen-subspace  and  the  cross-spectral  subspace.  The  resulting  plots  of 
Pd  versus  the  normalized  output  SINR  are  presented  in  Figure  12.  The  normalized  SINR  is 
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Figure  12;  The  AMF  probability  of  detection  versus  normalized  output  SINK. 

calculated  by  dividing  the  output  SINK  of  each  reduced-rank  detector  by  the  output  SINR  of 
the  full-rank  detector.  It  is  seen  that  the  full-rank  AMF  detector  and  the  rank  50  cross-spectral 
AMF  detector  demonstrate  very  similar  performance.  By  comparison,  the  rank  50  eigen-subspace 
detector  performs  quite  poorly.  Recall  that  a  target  input  SNR  of  15  dB  is  used  in  the  graphical 
examples  presented  earlier  in  this  section.  The  full-rank  and  cross-spectral  AMF  detectors  both 
yield  a  detection  probability  nearly  equal  to  unity  for  this  target.  The  eigen-subspace  detector, 
however,  yields  a  Pd  near  zero  for  this  same  scenario.  The  rank  50  cross-spectral  subspace  detector 
yields  a  12  dB  improvement  in  detection  sensitivity  over  the  rank  50  eigen-subspace  detector.  This 
corresponds  to  the  estimated  12  dB  gain  of  the  cross-spectral  subspace  processor  over  the  SINR 
obtainable  by  the  eigen-subspace  processor  at  rank  50  as  it  is  depicted  in  Figure  11. 
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8  Conclusions 


In  this  report  a  general  methodology  is  introduced  in  order  to  describe  reduced-rank  detectors 
and  to  analyze  their  statistical  properties  and  performance.  Although  the  AMF  CFAR  test  is 
emphasized  in  this  report,  the  analysis  technique  developed  herein  may  be  extended  readily  to 
nearly  any  other  test  by  simply  relating  the  variables  which  compose  the  test  to  the  GSC.  For 
example,  the  generalized  likelihood  ratio  test  (GLRT)  [23],  now  denoted  A3,  results  in  another 
CFAR  test  which  takes  the  form 


1  +  x^R  ^x 


A3  = 


=  1  + 


\yf 


1  +  xHR->x  -  ‘  "  P  (1  +  b«Rj->b) 

For  this  case  the  reduced-rank  GLRT  is  expressed  in  the  form, 


<  V3- 


A3,.  —  1  -b 


p  <  ^3) 

Ho 


(109) 


(110) 


where  the  term  Vr  =  ■s/pVr  and  p  is  the  reduced-rank  RMB  loss  factor  in  (96). 

The  computational  complexity  requirements  of  the  eigendecomposition  performed  for  subspace 
selection  suggest  that  other  procedures  need  to  be  discussed.  It  is  noted  that  other  unitary  operators 
U  may  be  used  with  the  cross-spectral  metric  to  determine  the  transformation  to  yield  new 
reduced-rank  tests  which  may  approximately  exhibit  the  optimal  performance  obtained  herein.  The 
Gram-Schmidt  procedure  is  one  example  of  a  lower-complexity  transformation  which  can  result  in 
excellent  performance  with  the  cross-spectral  metric.  Another  example  is  the  approximation  of  the 
cross-spectral  metric  which  uses  a  unitary  operator  that  does  not  completely  decorrelate  the  noise 
subspace  covariance  matrix.  This  technique  is  used  in  [19],  where  the  DCT  matrix  was  used  in 
order  to  reduce  the  computational  complexity. 

Finally,  it  is  possible  to  form  new  detectors  by  incorporating  other  constraints  in  conjuction 
with  the  STAP  detection  linear  constraint.  Examples  of  this  include  derivative  constraints  [34]  and 
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other  additional  linear  constraints  which  take  advantage  of  any  available  a  priori  knowledge  of  the 
interference  environment  [35,  36].  Fixed  transformations  of  the  data  prior  to  adaptive  filtering  also 
belong  to  this  category  and  the  popular  element-space  pre-Doppler  and  post-Doppler  algorithms, 
as  well  as  their  beamspace  analogies,  may  be  formulated  in  the  context  provided  by  this  report. 
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Appendix:  Algorithms  for  Finding  an  Orthonormal  Signal  Block¬ 
ing  Matrix  B 


In  this  appendix,  two  simple  algorithms  are  described  to  determine  an  orthonormal  signal  blocking 
matrix  B  given  the  normalized  N  X  1  dimensional  steering  vector  s.  The  first  algorithm  is  more 
accurate  and  utilizes  the  singular  value  decomposition  (SVD).  The  second  algorithm  is  faster,  and 
uses  the  QR  decomposition.  Both  of  these  algorithms  are  derived  and  explained  in  detail  in  many 
linear  and  numerical  algebra  texts,  such  as  [37]. 

Algorithm  1 

[U,S,,V]  =  svd(s'^') 

B  =  V(:,2:iV) 


Algorithm  2 


[Q,R]  =  qr(s) 

B  =  Q^(:,2:iV) 


In  many  cases,  a  non-orthonormal  signal  blocking  matrix  may  be  easily  constructed  and  the 
Gram-Schmidt  proceedure  can  be  implemented  to  provide  the  orthonormalization.  One  example 
of  this  is  a  system  which  utilizes  a  primary  space-time  steering  vector  to  set  the  array  manifold  in 
angle  and  Doppler.  The  normalized  steering  vector  s  is  then  fixed  to  point  broadside  with  respect 
to  the  array.  The  non-orthonormal  signal  blocking  matrix,  under  these  conditions,  is  implemented 
by  simply  taking  the  difference  of  the  appropriate  input  signals.  This  is  identical  to  the  array 
configuration  presented  in  [2]  for  narrowband  arrays. 
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MISSION 

OF 

ROME  LABORATORY 


Mission.  The  mission  of  Rome  Laboratory  is  to  advance  the  science  and 
technologies  of  command,  control,  communications  and  intelligence  and  to 
transition  them  into  systems  to  meet  customer  needs.  To  achieve  this, 
Rome  Lab; 


a.  Conducts  vigorous  research,  development  and  test  programs  in  all 
applicable  technologies; 

b.  Transitions  technology  to  current  and  future  systems  to  improve 
operational  capability,  readiness,  and  supportability; 

c.  Provides  a  full  range  of  technical  support  to  Air  Force  Materiel 
Command  product  centers  and  other  Air  Force  organizations; 

d.  Promotes  transfer  of  technology  to  the  private  sector; 

e.  Maintains  leading  edge  technological  expertise  in  the  areas  of 
surveillance,  communications,  command  and  control,  intelligence,  reliability 
science,  electro-magnetic  technology,  photonics,  signal  processing,  and 
computational  science. 


The  thrust  areas  of  technical  competence  include:  Surveillance, 
Communications,  Command  and  Control,  Intelligence,  Signal  Processing, 
Computer  Science  and  Technology,  Electromagnetic  Technology, 
Photonics  and  Reliability  Sciences. 


